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We have measured the proton and deuteron spin structure functions g’f and gf in the region of the 
2 2 12 12 12 

nucleon resonances for W <5 GeV and Q ~ 0.5 and Q ~ 1.2 GeV by inelastically scattering 
9.7 GeV polarized electrons off polarized ®®NH 3 and ®®ND 3 targets. We observe significant structure 
in g^ in the resonance region. We have used the present results, together with the deep-inelastic 
data at higher VF®, to extract r(Q®) = gi(x, Q^)dx. This is the first information on the low-Q® 
evolution of P toward the Gerasimov-Drell-Hearn limit at Q® = 0. 
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The nucleon spin structure functions gi and 52 have 
been and continue to be intensively studied via deep- 
inelastic lepton scattering Hi; however, they remain 
largely undetermined in the resonance region, where the 
squared invariant mass of the final state is less than 
4 GeV^. Here, gi and 52 can be used to investigate the 
helicity structure of the resonance transition amplitudes. 
States of definite spin and parity are more easily under¬ 
stood in terms of the virtual photon asymmetries Q| 


Ai{x,Q^) 

A2{x,Q^) 


g'1/2 ~ 0'3/2 _ J_ 
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in which is the squared four-momentum transfer, v 
is the energy transfer, x = /2Mv, M is the nucleon 

mass, and the structure functions Fi, gi and g 2 depend 
on both x (or W'^) and Q^. The cross sections ai /2 and 
cr 3/2 measure the strength of virtual transverse photon 
absorption leading to final-state spin projections of ^ and 
The cross sections ax = {o'i /2 + o' 3 / 2 )/ 2 , and 
measure longitudinal, transverse and interference photon 
absorptions. Positivity limits require that |Hi| < 1 and 
IH 2 I < y^R{x, Q2), in which R = ax/crT- The excitation 
of the A(1232) resonance (spin-|) allows for both i and | 
spin projections, and at low is expected to be primar¬ 
ily a magnetic dipole transition, for which a 2 , 12!<^112 = 3 
and Ai = —i. On the other hand, the 5'ii(1535) reso¬ 
nance has no spin-| projection, so Ai should be unity. 
The observed values of Ai in the resonance region are 
a combination of asymmetries for individual resonances 
and for the nonresonant background. The goals of the 
present measurements are to gain a better understanding 
of the resonances and to determine their influence on the 
deep-inelastic results, both in terms of radiative correc¬ 
tions and the evolution of r((3^) = /g gi{x,Q^)dx. Al¬ 
though the resonant contribution to P is insignificant at 
high it dominates the integral below about = 0.5 
GeV^. In fact, the limits set by the Gerasimov-Drell- 
Hearn (GDH) sum rule indicate that r((5^) should 
change sign in the region 0 < < 1 GeV^ and approach 

zero as ^ 0 . 

The data for the present analysis are part of the E143 
data set ||,P taken with a 9.7 GeV polarized electron 
beam (average polarization 85%) and cryogenic ^®NH 3 
and ^®ND 3 targets (average polarizations of 65% and 
25%, respectively). The data taken with the two spec¬ 
trometers at 4.5° and 7° corresponded to ~ 0.5 and 
1.2 GeV^ in the resonance region. 

Since these data were taken with longitudinal target 
polarization only, the determination of gi or Ai requires 
additional assumptions about either g 2 or A 2 . We ex¬ 
tract gi rather than Ai because gi is signihcantly less 
affected by our lack of knowledge of g 2 or A 2 . We set 
A 2 = 0 (corresponding to gi = —( 72 ) in our analy¬ 


sis. This is motivated by the fact that IA 2 I < '/R and 
existing data [|j indicate that R is small in the reso¬ 
nance region (R = 0.06 ± 0.02 for 1 < < 8 GeV^ 

and IT^ < 3 GeV^). We have explored the sensitivity 
to A 2 by considering the alternate possibilities g 2 = 0 
and 52 = 9 ^^ = —91 + gi(x')jx'dx', the Wandzura- 
Wilczek [ p^ twist-two form. Maximum deviations in gi 
from the A 2 = 0 case define the systematic errors due to 
uncertainty in A 2 . Even if A 2 were as large as 0.3, the ex¬ 
tracted values of gi would shift by less than 0.014, which 
is small compared to the statistical errors on each point. 

We have extracted gi from the absolute cross section 
differences for electron and nucleon spins aligned along 
the beam axis either parallel (tt) or antiparallel (ti) to 
each other @] : 


1 rf^Ti(TT) 

w 

^{E + E' cos 9 + i 2 )gi - ^/^ElA-, 


^ dndE' - T? + (f)^i + (f) 
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in which ctm is the Mott scattering cross section, E {E') 
is the initial (final) electron energy, 9 is the laboratory 
scattering angle, v = E — E\ and Ei and F 2 are the 
unpolarized structure functions. 

This method requires good knowledge of spectrome¬ 
ter acceptances, the number density of polarizable pro¬ 
tons or deuterons in the target, and detector efficiencies. 
Alternatively, one could extract gi from the count rate 
asymmetry as in Ref. |^,^; however, the dilution factor 
(the fraction of scatterings coming from a polarizable nu¬ 
cleon in the target) is needed in this case, which is more 
difficult to obtain reliably in the resonance region. Nev¬ 
ertheless, when we tried both methods, we found that 
they agreed to within a fraction of the statistical errors 
on each point (typically better than 3%). 

As the first step to determine the absolute cross sec¬ 
tion differences, we calculated the raw count differences 
per incident charge N^/qF — iVH/gll, with each term 
corrected for dead-time in the trigger electronics. For 
each electron, W'^ was calculated using the momentum 
and scattering angle determined from tracking. The data 
were binned in normalized by the product of target 
and beam polarizations and corrected for absolute spec¬ 
trometer efficiency. Each detector’s efficiency was de¬ 
termined by making a strict cut to select good electron 
events without using one of the detectors and checking 
how often that detector registered the electron. The ab¬ 
solute spectrometer efficiency is the product of all of the 
individual detector efficiencies (no evidence for correla¬ 
tions between them was found) |^. 

Fully corrected cross section differences dAa/dE'dH. 
were obtained with the help of a Monte Carlo simulation. 
This simulation was used to normalize the raw data for 
target density and spectrometer acceptance. The nor¬ 
malized data were then corrected for radiative and res¬ 
olution effects by an additive term determined from the 
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Monte Carlo routine. Small corrections for polarized 
and in the target (and in the case of the deuteron) 
were applied as in Ref. ||]. 

The Monte Carlo code simulated all relevant aspects 
of the experiment and was able to predict total count 
rates and count rate differences based on a set of input 
tables of cross sections and asymmetries. The unpolar¬ 
ized cross sections were calculated from a fit by Stuart 
et al. . The asymmetry tables were calculated from a 
combination of resonant and non-resonant contributions. 
The non-resonant part came from a parameterization of 
all existing deep inelastic data (Fit III of Ref. Q|), which 
was extrapolated into the resonance region. The asym¬ 
metry of the A(1232) resonance was fixed at Ai = —0.5. 
The contribution to the asymmetry from all other reso¬ 
nances was approximated by two constant values (one for 
the region below = 2.6 GeV^ and the other for the 
region above = 3.2 GeV^) and a linear interpolation 
between them. This simple parameterization was suffi¬ 
cient to achieve a reasonable fit to the data for the pur¬ 
pose of radiative and smearing corrections. The quality 
of the Monte Carlo results was also tested by comparing 
the predicted and measured total count rates and (quasi- 
) elastic asymmetries. We found good agreement within 
the statistical uncertainty of the data. 

The full model was used to calculate cross sections and 
asymmetries both in lowest order of QED (Born approx¬ 
imation) and with full radiative corrections following the 
prescription by Shumeiko et al. |Q. We used the peak¬ 
ing approximation after convincing ourselves that the 
full code produced negligible differences. The radiated 
cross sections and asymmetries were tabulated as input 
for the Monte Carlo code, while the Born cross sections 
and asymmetries were compared with the Monte Carlo 
output to determine the normalization factor and the ad¬ 
ditive correction for the raw data. 

Fig. 1 shows gi obtained from cross section differences 
for proton and deuteron (per nucleon) measured with the 
two spectrometers as a function of W^. The full length of 
the error bars corresponds to the statistical and system¬ 
atic uncertainties added in quadrature. The cross bars 
indicate statistical errors alone, which dominate the to¬ 
tal errors. Plotted as triangles are the data of Baum et 
al. 0 taken at similar kinematics and converted to gi 
for comparison by assuming A 2 = 0. Within errors, the 
two measurements agree well. The solid lines show the 
Monte Carlo simulation. The dashed curve is a calcu¬ 
lation using the code AO by Burkert and Li [|^. This 
calculation does not adequately describe the data above 
= 3 GeV^, because the non-resonant background 
contains only the single tt Born term. Although AO de¬ 
scribes the resonance region at higher rather well, 
there is a significant discrepancy in the second resonance 
region (2 < < 3 GeV^) at the lower 

Extracting the total (resonant plus non-resonant) Ai 


at fixed from the relationship gi = (Ai -I- 

^A 2 )Fi requires several simplifying assumptions. We 
assume A 2 = 0, i? = 0.25/Q^ for the non-resonant cross 
section 0 , and R = 0 for the resonant cross section 
and obtain Fi from measured values of F 2 and R. 
In that case, all data for both proton and deuteron are 
consistent with Af' = — ^ with a relative uncertainty of 
40-100%. This is also the asymmetry that AO predicts. 
On the other hand, Ai in the region = 2 — 3 GeV^, 
which is dominated by the 5'ii and D 13 resonances, is 
surprisingly high. Here the proton data at 4.5° and 7° 
are consistent with Ai = 0.9 ± 0.2, which is close to 
the positivity limit and well above the AO prediction. 
This could indicate a significant contribution from A 2 , 
from a non-resonant asymmetry that is much larger than 
expected from deep-inelastic extrapolations, or from a 
stronger than expected contribution from for which 
Al = 1. The smaller corresponding asymmetry for the 
deuteron (Ai = 0.3 ±0.2) may arise from Fermi smearing 
or a genuine n-p difference. More detailed information on 
the asymmetries of individual resonances and the non¬ 
resonant background are expected once high-precision, 
semi-exclusive data from TJNAF become available. 

Fig. 2 shows the integrals r(Q^) for proton and neu¬ 
tron, evaluated at the average for the resonance re¬ 
gion (M^ < < 4 GeV^). We summed our resonance 

results directly (where does not vary much) and then 
added a contribution from smaller x (larger W^) at the 
same fixed taken from Fit III to the world’s deep- 
inelastic gi data |^. The neutron integrals were derived 
assuming a 5% D-state probability for the deuteron. The 
statistical errors assigned to Fit III at given values of 
X and corresponded to the kinematically closest E143 
data points at 9.7 or 16 GeV, which dominated the fit 
in this region. Systematic errors were calculated using 
the systematic uncertainties for the measured gi in the 
resonance region added linearly to the systematic errors 
for the deep-inelastic region, which are highly correlated 
with each other. Extrapolation errors for the region be¬ 
low the last measured datum at a; = 0.03 were taken to 
be as large as the value that Fit III yields for x < 0.03. 

Although several models for the evolution of r(Q'^) 
exist pq-poll , we show here onl^wo representative ones, 
together with the evolution 1^ of the world’s deep- 
inelastic data due to the changing coupling constant as. 
Although the GDH sum rule is strictly valid only at 
— 0 where r(Q^) vanishes, it can be used to predict 
the slope of r(Q^) for small Q^. The solid line at low 
in Fig. 2 shows F = in which k is the anoma¬ 

lous magnetic moment of either the proton or neutron. 
Burkert and Ioffe ||^ consider the contributions from the 
resonances using the code AO, and the nonresonant con¬ 
tributions using a simple higher-twist-type form fitted to 
the deep-inelastic data. Their model is constrained to 
fit both the GDH and the deep-inelastic limits, and it 
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describes the data quite well. Soffer and Teryaev as¬ 
sume that the integral over gi + g 2 varies smoothly from 
high where 52 ~ 0 down to = 0. Using their sim¬ 
ple prediction for this integral and subtracting the con¬ 
tribution from g 2 alone using the Burkhardt-Cottingham 
sum rule gives the dashed curve in Fig. 2b, which 
also agrees quite well with our data. 

The present spin structure function data in the region 
of the nucleon resonances allow us to determine the in¬ 
tegrals r(Q^) for the first time at below 2 GeV^. In 
contrast to the nearly flat behavior in the deep-inelastic 
region above = 2 GeV^, F varies rapidly below 
= 2 GeV^. Models that interpolate between the deep- 
inelastic and GDH limits describe the data quite well in 
this non-perturbative regime. 
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TABLE 1. Integrals r(Q^) of the structure functions gi 
for the proton (p), deuteron (d) and neutron (n). The mea¬ 
sured sum Pres for the resonance region (W^ < 4 GeV^) 
is listed separately from the totals Ftot, which includes the 
deep-inelastic region as given by fits to the world’s data. 


(GeV^) 


Fres (±stat.±syst.) 

rtot(±stat.±syst.) 

0.5 

P 

0.026±0.008±0.008 

0.049 ±0.008 ±0.013 

0.5 

d 

0.000±0.013±0.008 

0.003 ±0.013 ±0.010 

0.5 

n 

-0.026±0.028±0.020 

-0.043 ±0.029 ±0.025 

1.2 

p 

0.040±0.003±0.004 

0.100 ±0.005 ±0.012 

1.2 

d 

0.026±0.006±0.004 

0.043 ± 0.008 ± 0.007 

1.2 

n 

0.016±0.013±0.010 

-0.006 ±0.019 ±0.020 
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TABLE II. Fully corrected cross section differences 
d{a^^ — /dE'd£l = Act (nb/GeV-sr) and structure func¬ 

tions gi for the proton in the resonance region. Listed are 
the data for the 4.5° and 7.0° spectrometers. The values of 
and (GeV^) are given at bin centers. The additive 
correction A(t_r includes both radiative effects and resolution 
corrections. 




Q" 

Act 

Ao-fl 

gi 

4.5° 

1.31 

0.55 

-74±133 

-714 

-0.013±0.023±0.023 


1.69 

0.54 

-289±116 

-244 

-0.070±0.028±0.006 


2.06 

0.53 

627±106 

282 

0.197±0.033±0.039 


2.44 

0.52 

828±109 

229 

0.320±0.042±0.046 


2.81 

0.50 

245±108 

-120 

0.113±0.050±0.025 


3.19 

0.49 

120± 80 

-166 

0.064±0.043±0.034 


3.56 

0.48 

75± 78 

-59 

0.046±0.048±0.015 


3.94 

0.47 

-2± 79 

-29 

-0.001±0.055±0.010 


4.31 

0.46 

158± 68 

-7 

0.123±0.053±0.015 


4.69 

0.45 

124± 66 

-10 

0.106±0.057±0.010 

7° 

1.19 

1.28 

-4± 5 

-19 

-0.003±0.003±0.003 


1.56 

1.26 

-8±16 

-49 

-0.007±0.014±0.005 


1.94 

1.23 

43±12 

4 

0.043±0.012±0.015 


2.31 

1.20 

162±16 

45 

0.193±0.020±0.031 


2.69 

1.18 

121±16 

3 

0.165±0.022±0.016 


3.06 

1.15 

95±15 

-12 

0.147±0.024±0.019 


3.44 

1.12 

57±15 

1 

0.098±0.026±0.009 


3.81 

1.09 

70±15 

2 

0.134±0.029±0.012 


4.19 

1.07 

110±15 

5 

0.233±0.031±0.025 


4.56 

1.04 

83±14 

4 

0.193±0.033±0.020 


4.94 

1.01 

63±14 

3 

0.158±0.036±0.022 


TABLE III. Same as Table 2, but for the deuteron. 





Act 

Ao-fl 

5i 

4.5° 

1.31 

0.55 

-195±207 

-397 

-0.034±0.036±0.021 


1.69 

0.54 

-190±185 

-211 

-0.046±0.045±0.015 


2.06 

0.53 

401±163 

80 

0.126±0.051±0.016 


2.44 

0.52 

127±165 

43 

0.049±0.064±0.007 


2.81 

0.50 

-41±145 

-34 

-0.019±0.067±0.005 


3.19 

0.49 

166±128 

-32 

0.089±0.069±0.013 


3.56 

0.48 

36±120 

-17 

0.022±0.074±0.004 


3.94 

0.47 

-8±112 

-14 

-0.006±0.078±0.002 


4.31 

0.46 

-27T101 

-9 

-0.021±0.078±0.003 


4.69 

0.45 

67±97 

-9 

0.058±0.084±0.007 

7° 

1.19 

1.28 

79±24 

18 

0.053±0.016±0.010 


1.56 

1.26 

1±21 

-30 

0.001±0.018±0.007 


1.94 

1.23 

30±19 

1 

0.030±0.019±0.007 


2.31 

1.20 

26±22 

11 

0.031±0.026±0.005 


2.69 

1.18 

28±22 

0 

0.038±0.030±0.004 


3.06 

1.15 

72±22 

0 

0.112±0.034±0.013 


3.44 

1.12 

19±22 

3 

0.034±0.038±0.004 


3.81 

1.09 

27±21 

2 

0.053±0.041±0.006 


4.19 

1.07 

42±21 

2 

0.089±0.044±0.010 


4.56 

1.04 

34±21 

1 

0.080±0.048±0.008 


4.94 

1.01 

50±20 

1 

0.126±0.051±0.015 
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Fig. 1. Measurements of gi as a function of for the proton at (a) 4.5° and (b) 7°; and for the deuteron at (c) 
4.5° and (d) 7°. The present data (circles) are plotted together with the data of Baum et al. (triangles), our Monte 
Carlo simulation (solid line), and the model AO of Burkert and Li (dashed line). The full error bars correspond 
to statistical and systematic errors added in quadrature, whereas the cross bars indicate statistical errors only. 


6 






































(GeV^) 


Fig. 2. Integrals of gi at several fixed values of for 
(a) the neutron and (b) the proton. The present data 
(circles) are plotted together with data from CERN 
(triangles), E143 deep-inelastic (squares), and E142 
1^ (inverted triangle). The curves correspond to the evo¬ 
lution Q of the deep-inelastic results due to changing 
as (solid line), the predictions of Burkert and Ioffe 0 ] 
(dotted line), the model of Soffer (dashed line), and 
the GDH approach to = 0 (solid line). Errors are 
indicated as in Fig. 1. 
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